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Abstract 

This paper deals with the multiobjective version of the optimal spanning tree problem. More precisely, we are in- 
terested in determining the optimal spanning tree according to an Ordered Weighted Average (OWA) of its objective 
values. We first show that the problem is weakly NP-hard. In the case where the weights of the OWA are strictly 
decreasing, we then propose a mixed integer programming formulation, and provide dedicated optimality conditions 
yielding an important reduction of the size of the program. Next, we present two bounds that can be used to prune 
subspaces of solutions either in a shaving phase or in a branch and bound procedure. The validity of these bounds does 
not depend on specific properties of the weights (apart from non-negativity). All these exact resolution algorithms are 
compared on the basis of numerical experiments, according to their respective validity scopes. 

Key words: Multiobjective spanning tree problem, ordered weighted average, MIP formulation, optimality 
conditions, branch and bound 



1. Introduction 

Multiobjective combinatorial optimization deals with problems involving multiple viewpoints [4j- More formally, 
the valuation structure of such problems is made of vectors (each component representing a specific viewpoint) instead 
of scalars. All popular single objective optimization problems (e.g., valued graph problems, integer linear program- 
ming...) can be recasted in this setting, and solution algorithms must be proposed. Two types of approaches can be 
studied: either one looks for a best compromise solution according to a given aggregation function (e.g., max operator, 



Chebyshev's norm to a reference vector 14011 . ordered weighted average [41], Choquet integral 11211 ). or one aims at 
generating the whole set of Pareto-optimal solutions (i.e., solutions that cannot be improved on one objective without 
being depreciated on another one). Computing this set (also called the Pareto set) is natural when no preferential 
information is available or when the available information is unsufficient to elicit the parameters of the aggregation 
function: the solutions in the Pareto set are indeed the only ones likely to be selected by any rational decision maker. 
The interest in this approach has spawned a substantial literature (for a survey on the topic, the reader can refer to sev- 
eral quite recent papersj5l0]). However, the number of Pareto-optimal solutions can grow exponentially with the size 



of the instance [e.g.[14j,l!5|] and the number of objectives 113 1H . The examination of the whole set of Pareto-optimal so- 
lutions can therefore become quickly awkward. Furthermore, focusing on one particular compromise solution makes 
it possible to considerably speed up the resolution procedure. Consequently, for practicality and efficiency reasons, 
the search for a best compromise solution seems to be a better approach when an aggregation function is available. 
This is the approach we study in this paper. 

More precisely, we investigate here the multiobjective version of the optimal spanning tree problem. This problem 
arises naturally in various contexts. For example, consider a broadcasting network where the values of the edges 
represent bandwidths. Assuming that the bandwidth of a chain equals the minimum bandwidth over its edges, it is 
well-known that, in a maximal spanning tree, the bandwidth between two nodes is the maximum possible. When there 
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are several scenarios of traffic (impacting the values of the bandwidths) or several opinions of experts on the values of 



the bandwidths, the problem becomes multiofr 
whole Pareto set in the biobjective case O 
are more than two objectives 



ective. Previous works on the subject mainly deal with generating the 
23>H3]> or computing a min-max (regret) optimal solution when there 
14411 . To our knowledge, there is therefore no operational algorithmic tool 
for this problem when there are more than two objectives and when the min-max (regret) criterion is not really suitable. 
The present paper precisely aims at tackling this gap, by providing algorithms able to optimize a less conservative 
decision criterion for any number of objectives. Provided the required preferential information is available, and 
provided the objectives are commensurate (which is the case in the above example for instance), we propose to resort 
to an averaging operator to compare the vectorial values of the feasible solutions (spanning trees). According to the 
decision context, one may however want to put the emphasis on the best, the worst or the median evaluations of a 
solution. In other words, one needs to assign importance weights not to specific objectives (scenarios, experts), but 
rather to best and worst evaluations. The Ordered Weighted Average (OWA) precisely enables to model such concern. 
For this reason, we focus here on the OWA-optimal spanning tree problem, i.e. finding the optimal spanning tree 
according to OWA in a multiobjective spanning tree problem. In the case of strictly decreasing weights (favou ring 
well-balanced solutions), the use of an OWA objective function has been studied by Ogryczak and Sliwinski l24l 12511 
in continuous optimization under linear constraints, and by Galand and Spanjaard ll ill in multiobjective heuristic 
search. In this paper, we propose new algorithms specifically dedicated to the optimization of OWA in multiobjective 
spanning tree problems, some of which are able to handle not only strictly decreasing weights but also every other 
types of weights. 

The paper is organized as follows. After recalling some preliminary definitions and stating the problem, we give 
some insights into computational complexity (Section 2). Then, we provide a mixed integer programming formulation 
of the problem, as well as a preprocessing procedure based on optimality conditions for the OWA-optimal spanning 
tree problem (Section 3). We then propose two - efficiently computable - alternative bounds for discarding sub- 
spaces of solutions in either a branch and bound or a shaving procedure (Section 4). Finally, we provide numerical 
experiments to assess the operationality of the proposed methods (Section 5). 



2. Preliminaries 



2.1. Multiobjective compromise search problem 

A multiobjective compromise search problem is a problem endowed with vectorial costs where one searches for 
a best compromise solution according to a given aggregation function. A generic multiobjective compromise search 
problem can be formulated as a mathematical program. We now introduce some notations for this purpose. We 
denote by X c {0, 1 } m the set of feasible solutions, / : X — > M. p a vector valued function on X, and (p : — > R a 
multiobjective aggregation function. Within this setting, a multiobjective compromise search problem is written as 
follows: 

I min <p(y) 
(P) s.ty = f(x) 
[ xeX 

Denoting by Y — f(X) = {/(x) : x e X) the image set of X in the objective space MP, problem P can be 
simply reformulated as min ve y <p(y). The OWA-optimal spanning tree problem obviously belongs to this class of 
problems. The resolution methods we propose hereafter for this problem are actually quite generic. The mathematical 
formulations of P will therefore be convenient to describe these methods in the sequel. We now more specifically 
introduce the OWA operator and then the OWA-optimal spanning tree problem. 



2.2. The OWA operator 

Given a set {1, . . . ,p) of objectives (to minimize), one can associate a vector in N p to every feasible solution 
of a multiobjective problem. The comparison of solutions amounts then to comparing the corresponding vectors. 
Following several works in multiobjective optimization [e.g. |23|, |27|, |28[], we propose to compare the vectors on the 
basis of their OWA value 14 111 (to minimize), defined as follows: 
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Definition 1. Given a vector y in W, its ordered weighted average is OWA(y)=Y? =i wiy®, where 2f =1 W; = 1 andy^i) 
> .. > y( p ) are the components of y sorted in non-increasing order. 

According to the decisional context, one can distinguish within the class of OWA operators two interesting sub- 
classes depending on the definition of the weights: 

• strictly decreasing weights, i.e. w\ > . . . > w p > 0: the set of weights naturally belongs to this class when 
performing robust discrete optimization. In robust optimization lfl9ll . the cost of a solution depend on different 
possible scenarios (states of the world). The aim is to find a robust solution according to this multiobjective 
representation, i.e. a solution that remains suitable whatever scenario finally occurs. The use of the OWA 
criterion in this setting is justified since it can be characterized by a set of axioms that are natural for modelling 



robustness 12711 . More precisely, these axioms characterize an OWA criterion with strictly positive and strictly 
decreasing weights. Compared to the max criterion frequently used in robustness, the OWA criterion is less 
conservative since it enables trade-offs between several scenarios. Note however that the OWA criterion includes 
the max criterion as a special case, when one sets wi = 1 - 2f =2 e ^ w 2 = £2, ■■ •> w p = s p with £2 > ■ ■ ■ > s p > 
and E2 tends towards 0. Another interesting special case is obtained for "big-stepped weights", i.e. when the 
gaps between successive weights are huge {w\ » . . . » w p ). The OWA criterion reduces then to the leximax 
operator, which consists in comparing two vectors on the basis of their greatest component, their second greatest 
one in case of equality on the first one, and so on... This criterion refines thus the max criterion by discriminating 
between vectors with the same value on the greatest component. 

• non-monotonic weights: one of the most famous decision criterion in decision under complete uncertainty 
(i.e., when several states of the world can occur, and no information is available about their plausibilities) is 
Hurwicz's criterion, that enables to model intermediate attitudes towards uncertainty (i.e., neither desperately 
pessimistic nor outrageously optimistic) by performing a linear combination of the maximum possible value of 
a solution under the different scenarios, and the minimum possible one. More formally, if y is the image of a 
solution in the objective space, then its value according to Hurwicz's criterion is: amaxiy, + (1 - a)min, ;y,. 
Hurwicz's criterion clearly is a special case of the OWA criterion, obtained by setting w\ = a, w p = 1 — a 
and Wj = Vi + I, p. As soon as there are more than two scenarios and a £ {0, 1}, these sets of weights are 
neither non-increasing nor non-decreasing. Another natural decision context where the sequence of weights is 
non-monotonic happens when every component represents the opinion of a particular expert, and one wants to 
dismiss the extreme opinions. For instance, one can set w\ = 0, w p = and w,- = l/(p - 2) V/ + I, p. 

2.3. Problem and complexity 

The problem we study in this paper is the OWA-optimal spanning tree problem, that can be formulated as follows: 

OWA-OPTIMAL SPANNING TREE PROBLEM (OWA-ST) 

Input: a finite connected graph G = (V,E), p integer valuations v e t for every edge e e E (2 = 1, . . . ,p), and a set of 
weights Wj{i - 1 , . . . , p) for the OWA criterion; 

Goal: we want to determine a spanning tree T* e arg min^eT OWA(/(T)), where T is the set of spanning trees in 
Gand/(r) = Q: (67 .vJ,...,2 <6 rvp. 

Coming back to generic formulation P of a multiobjective compromise search problem, problem OWA-ST cor- 
responds to the following specifications: a spanning tree T is characterized by a vector x = (jc 1 , . . . , x m ) e X where 
x-i = 1 when edge e' belongs to T, and its cost is f(x) = (f\(x), . . .,f p (x)) where /)(x) = Yj"L\ v\xK The aggregation 
function is of course defined by (p(y) = OWA(j). We now introduce a small instance of problem OWA-ST that will be 
used as a running example in the sequel of the paper. 

Example 1. Consider a clique G with 4 vertices, and assume that the edges have been evaluated according to 
3 objectives (scenarios, experts...). The set of vertices is V — {1,2,3,4) and the valuations of the edges are the 
following: v [U] = (3,2,3), v [1 ' 3] = (4,3,1), v [1 - 4] = (1,2,2), v [2 - 3] = (2,4,1), v [2 ' 4] = (2,6,1), v [3 - 4] = (1,5,1). 
A minimum spanning tree is T\ — {[1,4], [2, 3], [3,4]} according to the first dimension and the arithmetic mean, 
T2 — {[1,2], [1,3], [1,4]} according to the second dimension, and — {[1,3], [2,3], [3,4]} according to the third 
dimension. None of them is however completely satisfying: either it is too much unbalanced (f(T\) — (4,11,4) and 
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f(Ty) — (7, 12,3)), or it is too much conservative (/(TV) = (8,7,6)). Spanning tree T4 — {[1,2], [1,4], [2,3]} (with 
f(Ti\) — (6, 8, 6)) presents none of the drawbacks of the previous trees, and is OWA-optimal when setting for instance 
W\ — 0.5, W2 — 0.3 and W3 = 0.2. The overall ranking according to OWA on Ti,...,Tn is indeed: 1) T4 with 
CWA(6,8,6) = 7, 2) T 2 with CWA(8,7,6) = 7.3, 3) T x with OWA(4, 11,4) = 7.5, 4) T 3 with OWA(l, 12,3) = 8.7. 

As already indicated above, when w\ — 1 and all other weights w,'s are equal to zero, the OWA criterion reduces 
to the max criterion. In this case, it has been proved that computing a min-max spanning tree is NP-hard lfl4l l44ll . 
Consequently, problem OWA-ST is also NP-hard in the general case. Note however that it is polynomially solvable 
for some subclasses of instances: 

• when w\ = W2 = ■ ■ ■ = w p (arithmetic mean), an optimal solution can be obtained by valuing every edge e 
by Y?-=\ v ;> an d tnen a PPly m g a standard minimum spanning tree algorithm (e.g. Prim's algorithm or Kruskal's 
algorithm); 

• when w p = 1 and w, = Vi ^ p (min criterion), an optimal solution can be obtained by solving p standard 
minimum spanning tree problems (one for each objective) and then returning the optimal one among them 
according to OWA; 

• when there exists a permutation n of objectives such that v^ (1) > . . . > for every edge e, an optimal solution 
can be obtained by valuing every edge e by 2f =1 Wi v S,-)> an d then a Pply m g a standard minimum spanning tree 
algorithm. 

Note that these types of polynomial instances are quite uncommon. Nevertheless, one can design a Fully Polynomial 
Time Approximation Scheme (FPTAS) that works for any instance of OWA-ST. For this purpose, following Aissi et 
al. Jll], one can enrich the existing FPTAS for the determination of an approximate Pareto set in the multiobjective 
spanning tree problem l26Tl . This algorithm indeed returns a set T e of spanning trees, the cardinality of which is 
polynomial in the size of the instance. For every spanning tree T in the graph, there exists T e e T e such that 
f{T E ) < (1 + s)f(T) for i = I,..., p. Note that \y t < y\ V/] => OWA(y) < OWA(y') H, and that OWA((l + s)y) 
— (1 + e)OWA(y). Thus, a spanning tree T* which is OWA-optimal among the trees of T £ satisfies OWA(/(T*)) < 
(1 + e)OWA(/(r*)), where T* is an OWA-optimal spanning tree in T. We have therefore an FPTAS for problem OWA- 
ST (the algorithm is indeed polynomial since the search for T* is performed in a set of polynomial cardinality). The 
impact of this result is however mainly theoretical, due to the high computational complexity of the polynomial time 
procedure proposed by Papadimitriou and Yannakakis to compute an approximate Pareto set ll26Tl . For this reason, we 
provide in the following more operational exact algorithms (the complexity of which is exponential in the worst case) 
for problem OWA-ST. The validity of the proposed algorithms depends on the subclass of weights used in the OWA 
criterion: strictly decreasing or arbitrary (including the case of non-monotonic weights). 



3. Resolution method for strictly decreasing weights 

The strict decreasingness of the weights implies the convexity of the OWA function. It is well-known that convex- 
ity of the objective function is a nice property in minimization problems. In particular, this property makes it possible 
to linearize objective function tp = OWA in mathematical programming formulation P of Section 2, so as to get a 
mixed integer programming formulation of problem OWA-ST. This is the topic of the following subsection. 



3.1. MIP formulation 

In the way of Yaman et al. 14311 for the robust spanning tree problem with interval data, we start from a compact 
mixed integer programming (MIP) formulation of the minimum spanning tree problem proposed by Magnanti and 
Wolsey 12011 . In this formulation the minimum spanning tree problem is considered as a special version of a flow 
problem. Every edge e = [i, j] is replaced by two opposite arcs (i, j) and (J, i). The set of such arcs is denoted by A 
in the sequel. Some vertex of the graph — say 1 — is selected as a source, and n — 1 units of flow are incoming into 
it (where V = {1, ■ • • , n}). Furthermore, 1 unit of flow is outgoing of every vertex i + 1. To each feasible flow in the 
directed graph, one can associate a connected partial graph by selecting edge e = [i, j] as soon as there is at least one 
unit of flow on (i, j) or (j, i). By imposing that the number of selected edges is « - 1, one obtains a spanning tree. Let 
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ifij denote the flow on arc (/, j), x e denote the boolean variable taking value 1 if ifij > or <pjf > (for e = [i, j]), and 
v e denote the scalar value of edge e. The corresponding MIP formulation of the minimum spanning tree problem is: 



(^mst) 



nin ^ v e x e 

eeE 

(i,j)eA (j,i)eA 



77-1 if I = 1, 

-1 VieV\(l) 



t p ij <(n-l)x e Ve = [/,;], 
<Pji<(n-l)x e Ve = [i,j], 



1, 



<p>0, x e e{0, 1} Vee£. 

We modify and enrich this MIP formulation in order to take into account the multiobjective aspect of OWA-ST, 
as well as the OWA aggregation function. The multiobjective aspect can be modelled by simply adding constraints 
yi = TieeE v f^ V/ e jl, . . . , p) into Pmst, where y, is the variable corresponding to the value of the spanning tree 
on the i th component. The objective function is then OWA(j) = £f =1 w,y(,), yielding a new program denoted by 
Pov/a- Nevertheless, in this formulation, this program is not linear. A nice linearization of the objective function has 
been proposed by Ogryczak and Sliwinski [251. The key idea is to use the Lorenz vector of y, the i th component of 
which is defined, in a minimization setting, by L,(y) = Y!j=i y(f)- This notion has been introduced in economics for 
inequality comparisons (in a maximization setting), where y is seen as an income distribution [e.g. By noting 
that OWA(v) = Tr=i( w i ~ Wi + \)Lj(y), where w p+ \ = 0, it appears that function OWA is linear in the components of the 
Lorenz vector. Furthermore, component L,(y) is the solution of linear program P] defined below. However, when y 
is a variable, program P } is no more linear. To overcome this difficulty, it is worth considering the dual version D y ., 
where d'. (resp. r,) are the dual variables for the inequality (resp. equality) constraints: 



(P)) 



max 
s.t. 



0<o^<l V;e(l, 



,p}. 



mm 
s.t. 



If- , d l . 



n+d'j>yj V;e[l,...,H 
d\>0 Vye {!,..., p). 



Replacing Lfy) by program D y . for i = 1, 



(^owa) 



s.t. 



..,/?, program Powa becomes: 

p p 

'=1 7=1 

n + d'j^yj Vi,;e {l,...,p}, 

yi = IeeEVtx e Vie{l,...,p}, 

2 </>ij - Z ^ = |-i 

(i,j)eA {jfieA v 



if i = 1, 

V/€ V\{1] 



ifij < (n 
ifji < (n 



\)x e 
\)x e 



Ve 
Ve 



UJ], 



eeE 

d 1 > V/, j e { 1 , . . . , p) and r, unrestricted, 
<p>0, x e e{0, 1} VeeE. 

Note that this formulation is valid only when w,- - w, + i > for ; = 1, . . . ,p, which is the case here since the Wi's 
are assumed to be strictly decreasing. With m edges in the graph, the program involves p 2 + n + 2m + 1 constraints and 
p 2 + p + 3m variables (variables y,'s can be omitted in the implementation): its size is therefore linear in the size of 
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the input for a fixed number of objectives. In order to further reduce the number of variables and constraints involved 
in Pqwa, we now provide a preprocessing method that enables to reduce the density of the graph. 



3.2. Preprocessing of the graph 

As often done when presenting algorithms for constructing spanning trees [e.g. [37]], we describe our preprocessing 
phase as an edge coloring process. Initially all edges are uncolored. We color one edge at a time either blue (manda- 
tory) or red (forbidden). At each step of the preprocessing phase, we have therefore a partial coloring c(-) of the edges 
of G. For a partial coloring c, we denote by T{c) the set of trees including all blue edges, some of the uncolored ones 
(possibly none), and none of the red ones. The aim of the preprocessing phase is to color as many edges as possible, 
so as to get a maximal coloring c such that min7-<=7-( f ) OWA(/(T)) = min7- e 7- OWA(/(T)). We now give conditions 
under which an edge can be colored blue or red, which are adaptations of the well-known cut optimality condition and 
cycle optimality condition to our problem. 

Before introducing the optimality conditions, we recall some definitions from graph theory. A cut in a graph is a 
partition of the vertices into two disjoint sets and a crossing edge (with respect to a cut) is one edge that connects a 
vertex in one set to a vertex in the other one. When there is no ambiguity, the term cut will also be used to refer to the 
set of crossing edges defined by the partition of the vertices. 

Proposition 1 (optimality conditions). Let G be a connected graph with coloring c of the edges. The following 
properties hold: 

(i) Cut optimality condition. Let us consider a cut C in G with no blue edge. If there is some uncolored edge e e C 
such that OWA(v e — v e ) < Ofor all uncolored edges e' e C, then e belongs to an OWA-optimal tree in T(c). 
(//) Cycle optimality condition. Let us consider a cycle C in G containing no red edge. If there is some uncolored 
edge e € C such that OWA(v e — v e ) < for all uncolored edges e' € C, then e can be colored red without changing 
value mm TeT ( C ) OWA(f(T)). 

Proof. The proof relies on the following property of function OWA with non-increasing weights: OWA(y - y') > 
OWA(y) - OWA(y'). To prove this property, one uses the convexity of function OWA, which follows from the fact 
that OWA(y) = max^n 2; W0V(i) where II is the set of all possible permutations of (1, . . . ,p). Combining convexity 
of OWA and equality OWA(iy) = /tOWA(y), one deduces thus the required inequality: OWA(y - y') + OWA(y') = 
2(±OWA(y - y' ) + ±OWA(y')) > 20WA(±(y - /) + \y') = OWA(y). 

Proof of (J). Suppose there exists a cut C and an uncolored crossing edge e that satisfies the cut optimality 
condition. Let T be an OWA-optimal spanning tree of 7~(c) that does not contain e. Now consider the graph formed 
by adding e to T. This graph has a cycle that contains e, and this cycle must contain at least one other uncolored 
crossing edge — say e' — such that e' E C, and therefore we have OWA(v e - v e ') < 0. We can get a new spanning 
tree T e T(c) by deleting e' from T and adding e. We claim that OWA(f(T')) < OWA(/(r)). Cost f(T') is 
indeed equal to f(T) - v e + v e . By the property indicated above, we have OWA(/(T) - v e + v e ) - OWA(/(T)) < 
OWA(/(r) - v e ' +v e - f{T)) = OWA(v e - v e '). Consequently, since OWA(v e - v e ') < 0, we deduce OWA(/(r')) = 
OWA(/(r) - v e ' + v e ) < OWA(/(r)), and T is therefore an OWA-optimal spanning tree in T(c) (that does contain e). 

Proof of {fx). Suppose there exists a cycle C containing no red edge with an uncolored edge e e C such that 
OWA(/ - v e ) < for all uncolored edges e' e C. Let T be an OWA-optimal spanning tree of T(c) that contains 
e. Now consider the graph formed by removing e from T. This graph is compounded of two connected components. 
The induced cut contains at least one other uncolored crossing edge — say e' — such that e' e C, and therefore 
OWA(v e - v e ) < 0. We can get a new spanning tree T e T(c) by deleting e from T and adding e' . We claim that 
OWA(/(r')) < OWA(/(r)). Cost f(T') is indeed equal to f(T) - v e + v e ' . By the property indicated above, we 
have OWA(/(r) - v e + v e ') - OWA(f(T)) < OWA(v c ' - v e ). Consequently, since OWA(v e ' - v e ) < 0, we deduce 
OWA(/(r')) < OWA(/(T)), and T is therefore an OWA-optimal spanning tree in TTc) (that does not contain e). □ 

Let us come back to the instance of Example[TJto illustrate how to color edges according to these conditions. 

Example 2. In the clique of Example^ by considering cut {[1,4], [2,4], [3,4]}], edge [1,4] can be colored blue since 
CWA(v [M1 - v [2 ' 41 ) = CWA(-l,-4, 1) = -0.6 (with wi = 0.5, w 2 = 0.3 and w 3 = 0.2) and OWA(v [1A] - v [3 - 4] ) = 
OWA(0, -3, 1) = -0.1. Furthermore, by considering cycle {[2,4], [2, 3], [3,4]}, edge [2,4] can be colored red since 
CWA(v [231 - v [2 - 41 ) = OWA(0, -2, 0) = -0.4 and OWA(v [3A] - v [2 ' 4] ) = OWA{-\, -1,0) = -0.5. 
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Algorithm Preprocessing(G,c) 

Input : A graph G with coloring c of edges 

Output : A maximal coloring c 

for each edge e of E do 

if the cut optimality condition holds for e then 

set c(e) = blue 
else if the cycle optimality condition holds for e then 
set c(e) — red 
return c 



Figure 1 : Preprocessing algorithm in the case of strictly decreasing weights. 



Note that Pareto-dominance of edge e' over edge e implies OWA(v e - v e ) < since all components of v e ' - v e are 
negative in such a case (see edge [2, 4] in the previous example). Our optimality conditions are thus an enrichment of 
the multiobjective optimality conditions of Sourd and Spanjaard l35ll . which are sufficient conditions for an edge to 
be mandatory or forbidden when generating the whole Pareto set. These conditions are indeed exclusively based on 
Pareto dominance between edges. We take here advantage of the knowledge of the OWA operator to optimize so as 
to be able to color a higher number of edges. 



The single objective versions of these conditions make it possible to design a generic greedy method [e.g. I37H . 
from which Kruskal's and Prim's algorithms can be derived. For problem OWA-ST, this method can be adapted so 
as to preprocess the graph prior to the resolution of the MIP formulation by a solver. The corresponding algorithm is 
indicated in FigureQ] Obviously, and contrarily to the single objective case, this preprocessing method does not yield 
a spanning tree since binary relation OWA(v c - v e ) < does not induce a complete ordering of the edges in E. It 
enables however an important reduction of the number of variables (see numerical experiments), and therefore of the 
resolution time. 

The complexity of the preprocessing method strongly depends on the complexity of detecting uncolored edges 
satisfying an optimality condition. In practice, to determine whether an uncolored edge e = [i, j] satisfies the cut 
optimality condition, one performs a depth first search from i in the partial graph G£ ut = (V, Ef 1 ) where £^ llt = \e' e 
E | OWA(v e - v ) > 0} U [e' e E \ c(e') = blue}. If j does not belong to the set of visited vertices, then the partition 
between visited and non-visited vertices constitutes a cut for which e satisfies the cut optimality condition. Similarly, 
to determine whether an uncolored edge e = [i, j] satisfies the cycle optimality condition, one performs a depth first 
search from i in the partial graph GT = (V, E C J C ) where ET = W e E | OWA(v e ' - v e ) < 0} \ {e} U \e' € E \ c(e') = 
blue). If j is visited, then the chain from i to j in the search tree, completed with [i, j], constitutes a cycle for which 
e satisfies the cycle optimality condition. Since the number of edges in the graph is m and the complexity of a depth 
first search is within Oim) in a connected graph, the complexity of the overall preprocessing is 0(m 2 ). 

4. Resolution methods for arbitrary weights 

The linearization of the objective function is not valid anymore when the weights are not strictly decreasing. 
Facing the difficulty induced by the non-linearity of the objective function, we propose to resort to a branch and bound 
procedure when the weights are not necessarily strictly decreasing. For such procedures, the efficient computation of 
bounds on the optimal value of problem OWA-ST is therefore worth investigating. Besides, even when the weights are 
strictly decreasing, such bounds can be used in a shaving procedure (see Subsection 4.4 for a description) performed 
before the resolution phase in order to reduce the size of the instance. Thus, to obtain lower bounds, we present here 
two alternative relaxation methods for a multiobjective compromise search problem P: min^(y) s.t. y = fix), x e X, 
as defined in Subsectionl2.ll 



Relaxation of the image set: it consists in defining a subspace Y' 2 Y of the objective space (we recall that 
Y = f(X)) and solving problem P Y > defined by: 



(Py) 
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min ip(y) 
s.t. y e Y' 



• Relaxation of the objective function: it consists in defining a function <p' : MP —* M such that for all images y in 
Y, (f'(y) < ip(y) and solving problem P^ defined by: 



f min ip'(y) 
(P ¥ ) \ s.t.y = /(x) 
[ xeX 

The main point is then to define Y' (resp. ip') such that the optimal value of problem Py (resp. P^) can be 
efficiently computed and provide the tightest bound as possible. In this concern, it seems opportune to consider some 
continuous relaxations of Y (resp. linear aggregation function ip'). In the following, we show more precisely how 
these relaxations can be performed when the aggregation function is an ordered weighted average. Note that, in 
both relaxations presented in the following, the computational efficiency of the procedure (to solve Py or PS) only 
depends on the ability to quickly solve the single objective version of the problem. Consequently, these procedures 
can be applied to a broad spectrum of problems for which the single objective version can be solved efficiently. 

4.1. Defining Y' and solving problem Py 

We now present a first bound on the value of an optimal solution to P, obtained by relaxation of the image set. For 
this purpose, let us define relaxed space Y' in program Py as follows: 

Y' = {y e W : y* > b< Vi = 0,. ..,/>} 

where bj = min{f(x) : x e X) denotes the value of an optimal solution according to f, and fo = Yr=\ The values bfs 
(i = 1 , . . . , p) are obtained by n runs of a standard algorithm for the single objective version of the problem (provided it 
can be efficiently solved) for valuations fi's successively (i = 1 , . . . , p). Concerning the computation of bo, for simplic- 
ity purpose, we only detail the case of multiobjective spanning tree problems. The value bo is obtained by applying a 
standard minimal spanning tree algorithm on the graph where each edge e is valued by Vq = Yr=\ v " : f° r anv spanning 
tree characterized by x, we have indeed f Q (x) = £f =1 f(x) = Zf = i(Zee£ vfjt*) = Zeefi(I]f = i v e i )x e = 2 ee£ v e Q x e . Note 
that this technique to compute bo extends to many multiobjective versions of classical combinatorial optimization 
problems (actually, as soon as the value of a feasible solution is additively decomposable over its elements). 

As indicated above, in order for the relaxation to be interesting in practice, one needs to provide an efficient 
procedure to compute the optimal value of program Py. When the objective function is OWA, one clearly needs to 
take into account the way the components of solution y are ordered for solving Py . The OWA function is actually 
piecewise linear. More precisely, it is linear within each subspace of MP where all vectors are comonotonic, i.e. for 
any pair y,y' of vectors, there exists a permutation n of (1, . . . ,p) such that y„^) > . . . > y n ( P ) and > ■ ■ ■> y' n ( P y 
Problem Py> can thus be divided into several subproblems, each one focusing on a particular comonotonic subspace 
of W. The solution of Py reduces to solving each linear program Py<^ defined by a particular permutation n of 



(Py>,*) 



mrn£f =1 wovo 

s.t. y„ (i) > y„ (i+l) V;e{l,...,p-1}, (1.1) 

y t > bi V/e{l,...,p}, (1.2) 

2X1 J/ > bo, (1.3) 

y e W\ 



The value of the optimal solution y* Y , to Py is then min^n OWA(y Y , ) where y* Y , denotes the optimal solution to 
linear program Py^ and II the set of all possible permutations. Note that for p components, there are |II|=/?! linear 
programs to solve. However, in practice, it is not necessary to solve the pi linear programs. There exists indeed an 
easily computable permutation tt* for which OWA(y* Y ,) = OWA(y Y , »): 

Proposition 2 (Galand and Spanjaard, 2007 El). Let 71* denote the permutation such that b„-(i) > b„^2) > ■ ■ ■ > 
b„'( P ). For any feasible solution y to Py\ n , there exists a feasible solution y' to Pyx suc h mat OWA(y') — OWAiy). 

Proof. The idea is to determine a feasible solution y' to Py^- such that y..=y^ V/ (where y^\) > . . . > y( P ))- It im- 
plies indeed OWA(y) = OWA(y') and the conclusion is then straightforward. In this respect, we construct a sequence 



(y J ) j=i,...,k °f solutions and a sequence (ttO^i,,..,* of permutations such that yi is feasible for P Y , nl (for j=l, . . . ,k), with 
y'=y, n x -n, n k -n* and y'^y 2 ^ . . . -y k {j) V/. Assume there exist ;'o, i\ e {1, . . . ,p) such that i'o < i\ and b„\ {io) < b„\^y 
Let permutation n 2 be defined by tt 2 (Jo) = t 1 (ii)> T 2 (ii)=« 0°o)> and irii^n 1 ^) Vi + i'o, z'i. Let solution y 2 be defined by 
y 2 ,,. .=y ,,, for i=l, . . . , p. We now show that y 2 is a feasible solution to Py „2. Note first that y 1 .,. . > b n iti.\, y 1 ,,. . > 

^'(ii)' ^ii(i ) - ^ttH/,) anc * > ^t'('o)- Hence, constraints 1.2 are satisfied since: 
• ^(0=^.(0 ^ **»<0= ^ 2 (0 for ' * '0> »l- 

Constraints 1.1 are also satisfied since [y^ 1(i) > y' 1( . +1) V/] => [y 2 2(;) > y 2 , ( . +1) V/]. Indeed, we have ^, (l) =^2 (0 V/. 
These equalities imply also that constraint 1.3 is satisfied and that y 2 ;) = y^ V/. Solution y 2 is therefore feasible for 
P Y > y- with y 2 ^=y(i) V/. Since any permutation is the product of elementary permutations, one always can construct in 
this way a sequence of permutations that leads to n* (and the corresponding feasible solutions). By setting y'-y k , one 
obtains the desired feasible solution to Pyy ■ Q 

An immediate consequence of this result is that OWA(y^,) = OWA(yJ,, A Thus, the computation of OWA(y^,) 
reduces to solving linear program Pyx ■ For the sake of illustration, we present below an example. 

Example 3. Let us come back to Example \l\whe re bounds bi's are defined from spanning trees T\, T2 and T3 by 
b\ = 4, &2 = 7, bi — 3. The optimal spanning tree with respect to fo is spanning tree T — {[1,4], [2, 3], [3,4]} with 
fo(T) — 19. The value of bo is then 19. It yields the following program Py ( the values ofw\, W2 and w 3 do not matter): 

min OWA(y) = w\y m + w 2 y (2) + w 3 y (3) 
s.t. yj > 4 y 2 > 7 y 3 > 3, 

yi +y 2 + y 3 > 19, 

yi e R, y 2 e R, y 3 € R. 

As a consequence of Proposition^ solving Py amounts to solving linear program Py,n' defined by: 

min wiy 2 + w 2 yi + w 3 y 3 
s.t. y 2 >yi >y 3 , 

yi > 4 y 2 > 7 y 3 > 3, 

yj +y 2 +y 3 > 19, 

yi e R, y 2 e M, y 3 e R. 

We wave indeed n* (Y) — 2, tt*(2) = 1 and 7r*(3) = 3 since bi>b\> £> 3 . When w\ — 0.5, w 2 = 0.3 anc/ w 3 = 0.2, f/ze 
va/Me of the bound provided by programm Py,n' is therefore 6.5 (we recall that the value of the OWA-optimal spanning 
tree is 1 ). 

4.2. Defining tp' and solving problem P^i 

The second bound for P is obtained by relaxation of the objective function. The idea is to use a linear function ip' 
defined by <p' (y) = hyi where the vector A = (A\, . . . , A p ) of coefficients satisfy the following constraints: 

|/| 

/t,<^w, for all subset / c {1, p] of objectives (1) 

iel i=l 

When using such weight vectors, we have indeed <p'(y) < (p(y) for all y in Y, as established by the following: 

Proposition 3. Let A denote the set of all vectors A e W satisfying Constraint^ If A e A, then for all y e W, 
2f =1 A, yi < OWA(y). 

Proof. Let n denote a permutation such that y„^ > . . . > y„( P ). Note first that OWA(y) = Yr=\ woVo = 
2f = i(Zij'=i w j)(yn(i) - yn{i+\)) where y n {p+\) = 0. By Constraint[T]with subset of objectives / = {n(l), . . . , n(i)}, we have 
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£>i wj > £' j= i Mj) v/ - Thus Zf=i(E;=i w^)Cy„(i)-^(i + i)) > Zf=i(Z>i ^<o)Ov<)-3Vo+i))' Since 2f =1 (Z}=i -WOV*)- 

M'+d) = 2f = i ^(o^(o and 2f = i ^oMo = 2f=i we have owa(j) > £f =1 hyi- □ 



This result can be seen as a specific instanciation of a more general result of Schmeidler 11321 13311 on Choquet 
integrals in decision under uncertainty (OWA being a particular subclass of Choquet integrals). Schmeidler's result 
has been used to provide bounds in Choquet-based optimization under uncertainty Jg], as well as under multiple 
objectives J^]. Note that, in the case of decreasing weights in OWA, a useful by-product of Schmeidler's result is the 
existence of a normalized set of weights (i.e., summing up to 1) such that Proposition [3] holds: for instance, when 
setting Aj - l/p, we have ^i(l/p)yi < OWA(y) by Chebyshev's sum inequality. Nevertheless, when the weights are 
not decreasing, the existence of normalized weights satisfying Constraint[T]is not guaranteed. 

Solving problem P^ can be efficiently done by valuing every edge e by Yr=\ Mv\ ar, d men performing a standard 
minimum spanning tree algorithm. Besides, in order to obtain the best lower bound as possible, the optimal set 
of weights A e A (i.e., providing the best lower bound according to Proposition can be obtained by solving the 
following program: 

/' 

max z(A) = min ) A{yi, 

i=l 

|/| 

s.t. Yj Ai -Yi Wi w eu.-- ■>/>}. 

m i=i 

Ai>0 Vi=l,...,/>. 

Given that z is a concave piecewise linear function of A for a fixed y (since it is the lower envelope of a set of linear 
functions (2f =1 : y e I'D, we solve this program by using the SolvOpt library 1 17], which is an implementation 
of Shor's r-algorithm l34ll . This algorithm is indeed especially convenient for non-differentiable optimization, and 
the implemented SolvOpt library enables to perform constrained optimization. This approach is closed to the lower 
bounding procedure proposed by Punnen and Aneja l29ll for min-max combinatorial optimization. In broad outline, 
it can be viewed as a sequence of minimum spanning tree computations according to a varying weight vector A, until 
convergence towards a weight vector maximizing min ye Y Sf=i ^yi- 

Example 4. Let us come back to the instance of Example\l\ Running Shor's r-algorithm yields A — (0.28, 0.5, 0.22) 
and y — (6,8,6) (the image of spanning tree {[1,2], [1,4], [2,3]}, which is actually OWA-optimal). The bound is 
therefore 0.28 x 6 + 0.5 x 8 + 0.22 x 6 = 7 (which is the value of an OWA-optimal tree). This is better than the 
bound obtained by the previous relaxation (value 6.5). The computational burden is however more important, since it 
requires to solve much more single objective problems. 



4.3. Branch and bound procedure 

The two bounds presented above can of course be inserted into a branch and bound procedure in order to determine 
an OWA-optimal spanning tree for arbitrary weights. We summarize below the main features of the branch and bound 
procedure we propose. 

Branching scheme. The branching scheme is very simple: at each node, an edge e of G is selected and two subproblems 
are created. In the first one, edge e is mandatory while in the second one, edge e is forbidden. The heuristic to select 
the next edge consists in searching for an edge e such that Ya-=\ v 1 * s minimal among the remaining ones. 

Computing the lower bound. The lower bound at each node of the branching tree is computed by relaxation of the 
image set (by generating and solving problem Py>) or by relaxation of the objective function (by solving a sequence 
of problems P^). 

Updating the incumbent. When defining and solving problem Py> or P^, the algorithm checks whether a newly 
computed spanning tree is better than the current best known spanning tree according to OWA. We take indeed 
advantage of the property that feasible spanning trees are generated when computing the lower bound: either p + 1 
spanning trees are generated to obtain the values of bj (i = 0, . . . , p) in Py> or a sequence of spanning trees is generated 
during the running of Shor's r-algorithm to find the best possible bound according to P^. 
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Initialization. A branch and bound algorithm is notoriously more efficient when a good solution is known even before 
starting the search. The initial solution in our method is the best known solution (according to OWA) after the run of 
the shaving procedure described below. 

4.4. Shaving procedure 

The definition of lower bounds makes it possible to resort to a shaving procedure in order to reduce the size of 
the instance before running the main algorithm. The term "shaving" was introduced by Martin and Shmoys 12 ill for 
the job-shop scheduling problem. Assuming at least a feasible solution is known, this procedure works as follows: 
for each edge e, we build a subproblem in which e is made mandatory. If the computation of the lower bound proves 
that the subproblem cannot improve the current best known solution, then it means that e can be made definitively 
forbidden (colored red). Conversely when e has not been colored red by the previous procedure, we test similarly 
whether e can be made definitively mandatory (colored blue). Note that, here again, when computing the lower 
bounds, one checks whether a newly detected spanning tree improves the current best known solution (according to 
OWA). Of course, the shaving procedure is all the more efficient as a good feasible solution is initially known. For this 
purpose, we generate k feasible solutions by running a £-best ranking algorithm (i.e., returning the k best solutions) 
for the minimum spanning tree problem on the instance valued by the arithmetic mean of the vectors. The choice of 
k depends on the size of the instance. 

Example 5. Let us come back again to the clique of Example\l\ Assuming that the k-best ranking algorithm is run 
for k — 2, the current best known solution is then spanning tree T4 — {[1,2], [1,4], [2,3]}. Let us now simulate 
the progress of the shaving procedure. Making edge [1,2] forbidden, the minimum spanning tree for weight vector 
A — (0.3,0.5,0.2) is {[1,3], [1,4], [2,3]} with image y — (7,9,4). The induced lower bound is 7.4, which is strictly 
greater than the OWA-value of T4 (OWA{f(T^)) — 1). Therefore edge [1,2] is colored blue. Then making edge 
[1,3] mandatory, the minimum spanning tree for weight vector A — (0.5,0.3,0.2) is {[1,2], [1,3], [1,4]} with image 
y — (8, 7, 6). The induced lower bound is 7.3, which is strictly greater than the OWA-value 0/T4. Therefore edge [1,3] 
is colored red. Next, making edge [1,4] forbidden, the minimum spanning tree for weight vector A — (0.3,0.5,0.2) 
is {[1,2], [2,3], [3,4]} with image y — (6, 11,5). The induced lower bound is 8.3, and edge [1,4] is therefore col- 
ored blue. Finally, making edge [2,3] forbidden, the minimum spanning tree for weight vector A — (0.2, 0.5, 0.3) is 
{[1, 2], [1,4], [3,4]} with image y — (5,9, 6). The induced lower bound is 7.3, and consequently edge [2, 3] is colored 
blue. All edges ofT\ are colored blue: one concludes that T4 is an OWA-optimal spanning tree without even starting 
the main algorithm. 



5. Experimental results 

Before we study more carefully the behavior of our algorithms, we give some insights into previous results on 
related topics. The most widely studied related topic is the generation of the Pareto set in the bi-objective spanning 
tree problem To our knowledge, the most efficient algorithm for this problem enables to solve 

randomly drawn clique instances containing up to 400 vertices l35ll (note that these results significantly improved 
the size of the instances that could be handled, since it grew from 40 to 400 vertices). However these results do not 
extend to more than two objectives. More generally, even when looking for a single compromise solution within 
the Pareto set, it seems that there is very few available numerical experiments in the literature for more than two 
objectives. Although several works deal with the min-max spanning tree problem (i.e. determining a spanning tree 



minimizing the max criterion) |jj,|2j,[14j,[39j|44j], the content of these works is indeed mainly theoretical. Actually, the 
only numerical results we know for more than two objectives are devoted to the determination of a Choquet-optimal 
spanning tree 19J]. The size of the tackled instances goes from 30 to 70 vertices according to the number of objectives 
and the parameterization of the Choquet integral. 

5.1. Experimental details 

All the algorithms have been implemented in C++ and were run on a 2.6 GHz personal computer with a memory 
of 3.4GB. The test instances are defined as follows. All considered graphs are cliques. The components of cost vectors 
are randomly drawn between 1 and 100 on each edge. The number of objectives varies from 3 to 10, and the number 
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of vertices from 10 to 100 for strictly decreasing weights, and from 10 to 25 for non-monotonic weights. For each 
kind of instances (depending on the number of nodes and on the number of objectives), 30 instances were randomly 
drawn to obtain average results. 



5.2. Tests with strictly decreasing weights 

The global procedure to solve problem Powa when the weights are strictly decreasing consists of two phases: 

1. Coloration phase: making the most possible edges blue (mandatory) or red (forbidden) by running first the 
preprocessing procedure, and then a shaving procedure taking into account the coloration obtained by prepro- 
cessing. 

2. Resolution phase: determining the OWA-optimal spanning tree by running the main resolution algorithm (so- 
lution by MIP or branch and bound) on the reduced instance. 

In this section, one summarizes the results one has obtained by running this global procedure. For initializing the 



shaving procedure in the coloration phase, the £-best ranking algorithm proposed by Katoh et al. [18] is launched for 
k varying from 500 to 5000 depending on the size of the instance. In all cases, the order of magnitude of its running 
time is about a few milliseconds. Concerning the shaving itself, preliminary results have shown that it is much more 
efficient when the bound is defined by relaxation of the objective function than by the relaxation of the image set, 
especially when the size of the instances grows. For this reason, one only summarizes here the execution times 
obtained when the shaving procedure is performed by relaxation of the objective function. This shaving procedure 
is denoted by sh^ in the sequel. Finally, one compares the various resolution algorithms proposed in this paper: the 
solution by MIP and the one by branch and bound (two versions according to the bound adopted). The mixed integer 
program is solved by using solver ILOG CPLEX 1 1 . Regarding the branch and bound procedure, the bounds proposed 
in Section 4 are compared: the one obtained by relaxation of the objective function, and the one obtained by relaxation 
of the image set. 

Tables [T] |2] and |3] summarize the results obtained by running the algorithms on cliques with respectively 3, 5 or 
10 objectives (the weights of the OWA operator are indicated in the caption of the tables), for various numbers of 
vertices for which the average resolution time is lower than 30 minutes. The upper part of the tables summarizes 
the informations about the coloration phase (line pp for the preprocessing procedure and line s/y for the shaving 
procedure). For each procedure and each size of instance, the average execution time is indicated. Furthermore, 
below line pp (resp. shy/), the average number of edges made blue (#b) or red (#r) after preprocessing (resp. after 
preprocessing and shaving) is indicated (couple (#b - #r)). In the lower part of the tables are indicated the average 
total resolution times for various combinations of procedures for the coloration phase and the resolution phase. To 
evaluate the variability of the resolution time, the minimal running time (min) as well as the maximal one (max) 
are also indicated (couple min - max under the average execution time). To encode the combinations, the following 
abbreviations are used: MIP stands of course for Mixed Integer Programming, BBy for the branch and bound obtained 
by relaxation of the image set, and BB^ for the branch and bound obtained by relaxation of the objective function. 
Note that, in order to show the impact of the coloration phase, the first line of the lower part indicates the resolution 
times when solving the MIP formulation without resorting to any preprocessing or shaving. 

Table 1 shows that the best results are obtained by using the MIP formulation with a coloration phase. This is 
the only algorithm able to handle instances with more than 60 vertices with an average computation time below 30 
minutes. It is essential to note that the coloration phase has a very significative impact on the computation time. For 
n = 60, one sees indeed that the resolution time for MIP is above 1000 seconds, while it is below 100 seconds when 
the coloration phase is used. This can be easily understood by observing that the number of colored edges at the end 
of the coloration phase is a low fraction of the initial number of edges. For instance, in average, it remains only 67 
uncolored edges over 4950 for n = 100. The main interest of preprocessing is that it enables to color a lot of edges 
in a very low computation time. Hence, the number of edges that are tested during the shaving procedure is reduced, 
which is computationally interesting since the shaving is more powerful but takes also much more time. Performing 
shaving after preprocessing makes it possible to color blue a large part of the optimal spanning tree in a reasonable 
computation time. For illustration, in average, preprocessing colors red 4520.8 edges for n = 100, and shaving colors 
blue 71.6 edges (including those colored blue by preprocessing) over a maximum number of 99. A large amount of 
the total resolution time is spent in the coloration phase: for n = 100, 988.84 seconds are spent in average in the 
coloration phase (pp+sly) over a total resolution time of 1167.83 seconds. Finally, note that there is an important 
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n 


20 


30 


40 


50 


60 


PP 





0.02 


0.05 


0.11 


0.21 




(2.8 - 128.3) 


(4.4 - 332) 


(5.8 - 628.3) 


(6.3 - 1025.4) 


(7.1 - 1524.9) 


sly 


0.51 


2.89 


10.76 


31.43 


79.19 




(12- 160.7) 


(20.1 - 390) 


(26.1 -720.5) 


(33.5 - 1154.7) 


(37.9 - 1680.4) 


MIP 


0.51 


2.8 


8.11 


56.11 


1018.5 




0.02 - 4.65 


0.2 - 20.3 


0.45 - 49.7 


0.88 - 639.9 


1.58 - 16842 


pp+sly +MIP 


0.55 


3.04 


11.73 


32.26 


89.38 




0.31-0.99 


2.26 - 4.93 


7.88 - 32.55 


25.7 - 43.4 


68.8 - 177.3 


pp+sly +BB r 


0.58 


3.11 


25.89 


34.51 


241.52 




0.31 - 1.49 


2.26 - 5.79 


7.87 - 345.4 


25.7 - 67.7 


68.8 - 1945.5 


pp+sly +BB,/ 


1.35 


10.27 


72.1 


206 


1049.7 




0.32 - 3.79 


2.71-37.2 


8-571.2 


28.5 - 1183.3 


75.8 -4457.9 



n 


70 


80 


90 


100 


PP 


0.34 


0.54 


0.79 


1.14 




(8.5 - 2125.2) 


(9.3 - 2827.8) 


(12- 3624.1) 


(13.2-4520.8) 


sly 


165.6 


314.7 


563.7 


987.7 




(50.1 -2318) 


(52.8 - 3044.7) 


(63.1 - 3879.1) 


(71.6-4811.4) 


pp+sly +MIP 


266.94 


342.79 


590.98 


1167.83 




138.3 - 3206.5 


267.1 -720.6 


475.2 - 1025.3 


810.6-4453 



Table 1: Synthesis of the numerical results for 3 objectives (wi = 0.6, wj = 0.3 and »'3 = 0.1). Execution times are indicated in seconds. 



variability in the resolution time according to the instance (for n = 60, the MIP algorithm takes between 1.58 seconds 
and more than 4 hours and half). The coloration phase tends however to reduce this variability. 



n 


10 


20 


30 


40 


50 


60 


PP 





0.01 


0.02 


0.065 


0.14 


0.27 




(0.16-8.36) 


(0.18 -63.42) 


(0.27 - 187.2) 


(0.16- 384.9) 


(0.27 - 668) 


(0.3 - 1034.9) 


sly 


0.13 


1.71 


10.5 


43.57 


132.21 


337.13 




(3.54 - 26.84) 


(5.42 - 142.5) 


(8.13 -359.3) 


(9.92 -667.1) 


(12.77- 1079.7) 


(15.52- 1596.3) 


MIP 


0.07 


1.74 


52.39 


133.8 


1426.7 






0.03-0.14 


0.12 - 12.06 


0.59 - 1166.2 


0.94 - 3591 


12.26 - 16610 




pp+sry+MIP 


0.16 


2.63 


27.27 


87.66 


433.8 


1202.7 




0.09 - 0.24 


1.47 - 12.83 


8.66 - 284.2 


38.88 - 693 


110.5 -2589.6 


301.8 - 10153.4 


pp+sry+BB r 


0.17 


14.23 


764 










0.09 - 0.62 


1.5-218.8 


8.75 - 13767.5 








pp+sly +BB,/ 


0.37 


19.63 


541.3 










0.09 - 2.02 


2.58-91.58 


13.52-5018.5 









Table 2: Synthesis of the numerical results for 5 objectives (w\ = 0.5, W2 = 0.3, W3 = 0.1, w>4 = 0.06 and 11/5 = 0.04). Execution times are indicated 
in seconds. Symbol "-" means that the average execution time exceeds 30 minutes. 

Tables |2] and [3] show that the number of objectives has a great impact on the performances of the algorithms. For 
instance, for n - 60, the resolution time of algorithm pp+sty +MIP is 1202.7 seconds in average for 5 objectives, 
compared to 89.38 seconds for 3 objectives. Besides, for 10 objectives, one can only handle instances with up to 35 
vertices within the time limit of 30 minutes. The resolution time is also sensitive to the weights of the OWA operator. 
In order to evaluate this sensitivity, we also performed a test with weights inducing an OWA operator closer to the 
arithmetic mean. Namely, we set w\ — 0.4, wo = 0.35 and W3 = 0.25. 

For this set of weights, all the tested algorithms perform better, as it can be seen in Table [4] In particular, the 
preprocessing procedure makes it possible to color blue more edges than for the previous set of weights. For example, 
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n 


10 


15 


20 


30 


35 


PP 








0.01 


0.03 


0.05 




(0 - 2.22) 


(0 - 8.73) 


(0.03 - 20.2) 


(0.03 - 69.5) 


(0.03 - 108.5) 


sly 


2.69 


5.21 


11.41 


44 


86.57 




(3.78 -25.58) 


(2.43 - 58.33) 


(1.83 - 107.07) 


(1.1 -256.7) 


(1.3 - 365.27) 


MIP 


0.15 


1.87 


11.47 


519.2 






0.05 - 0.55 


0.13-7.65 


0.6 - 33.65 


4.76 - 6623 




pp+sty +MIP 


2.77 


6.54 


19.29 


311 


1284.8 




2.32-3.58 


4.42 - 12.52 


10.2 - 40.6 


41.26-2937.3 


188.1 -8504.6 


pp+sty +BBy< 


3.15 


91.74 










2.32 - 8.47 


4.41 -532.71 








pp+sfy +BB V , 


7.93 


58.27 


477.4 








2.32 - 25.75 


5.21 - 317 


13.32 -2109.5 







Table 3: Synthesis of the numerical results for 10 objectives (wi = 0.25, W2 = 0.2, w$ = 0.15, W4 = 0.1, W5 = 0.09, = 0.08, w-j = 0.06, 
W8 = 0.04, 11/9 = 0.02 and wio = 0.01). Execution times are indicated in seconds. Symbol "-" means that the average execution time exceeds 30 
minutes. 



n 


40 


50 


60 


70 


80 


PP 


0.06 


0.12 


0.22 


0.38 


0.59 




(23.57 - 717.6) 


(30.47- 1148.5) 


(35.9 - 1675.2) 


(41.85 - 2305.7) 


(48 - 3033.6) 


sty 


2.75 


7.22 


18.65 


37.51 


74.3 




(32.5 - 733.9) 


(41.9-1167.3) 


(49.3 - 1700.2) 


(57.1 -2334.9) 


(64.85 - 3064.8) 


MIP 


7 


65.24 


243.8 








0.25 - 94.75 


0.57 - 830.3 


7.12- 3227 






pp+sby +MIP 


2.81 


7.35 


18.9 


37.96 


75.21 




2.04 - 3.96 


4.22- 11.37 


14.23 - 22.75 


26.26 - 48.95 


55 -95.15 


pp+sry+BB F 


2.82 


7.38 


18.97 


38.376 


77.85 




2.04 - 4.01 


4.22- 11.86 


14.23 - 23.04 


26.26 - 55.59 


54.99- 104.1 


pp+sty +BB V ' 


5.37 


13.28 


51.5 


114.9 


391.1 




2.1 - 17.03 


4.33 - 48.46 


15.57 - 184.5 


26.57 -484.1 


64.08 - 1853.1 



Table 4: Synthesis of the numerical results for 3 objectives (ivi = 0.4, w% = 0.35 and W3 = 0.25). Execution times are indicated in seconds. 

for n = 60, 49.3 edges are colored blue in average (over a maximal number of blue edges of 59) while only 37.9 are 
colored blue in average when the previous set of weights is used. Moreover, one can also observed that, once the 
coloration phase is performed, branch and bound procedure BBj" is as efficient as the resolution by MIP. 

5.3. Tests with non-monotonic weights 

We present here the results obtained when the weights are non-monotonic. In this subsection, we used Hurwicz's 
criterion as OWA operator. We recall that, for a vector y, it is defined as amax, y, + (1 - ff)min, ;y,. The instances 
are defined similarly to the previous subsection, with 3, 4 or 5 objectives, and parameter a varying from 0.4 to 0.6 
in Hurwicz's criterion. As mentioned in Section 3, the MIP formulation is no more valid in this case, as well as 
the preprocessing procedure (the optimality conditions do not hold anymore). Furthermore, the bound obtained by 
relaxation of the objective function becomes weak since there does not necessarily exist a normalized set of weights 
satisfying Constraint Q] (see Section l4~2b . Consequently, in the global procedure described in Subsection 5.2, the 
coloration phase consists of a single run of the shaving procedure with the bound obtained by relaxation of the image 
set, and the resolution phase consists of applying branch and bound BBy. Tables [5] [6] and [7] summarize the results 
obtained. The conventions are the same than in the previous subsection. 

In these tables, one can observe that parameter a has a significant impact on the resolution times: it is easier for 
values of a that are above 0.5. For example, with 3 objectives and n = 25, the average execution time is 5.09 seconds 
for a = 0.6, while it is 195.6 seconds for a = 0.5, and it is more than 30 minutes for a = 0.4. The sizes of the 
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Yl 




1 


i j 








Shy/ 








0.02 


0.04 




a = 0.4 




(1.5 - 3.23) 


(0.13 - 12.53) 


(0-4) 


(0 - 2.73) 






Shy/ +BBy/ 


u 


U.Uj 


i n 1 
1.U1 


5 l.yi 








- 0.02 


- 0.77 


0.05-6.18 


0.85 - 129.7 






Shy/ 








0.02 


0.04 


0.07 


a = 0.5 




(1.72- 3.22) 


(0.07 - 10.83) 


(0.03 - 18.5) 


(0-15.27) 


(0 - 12.57) 




Shy/+BBy- 





0.03 


0.87 


7.53 


195.6 






- 0.03 


0-0.14 


0.06 - 9.93 


0.23 - 40.46 


2.14-2203.1 




Shy/ 





0.01 


0.01 


0.03 


0.06 


a = 0.6 




(1.73-3.3) 


(1.9-23.2) 


(0.8 - 48.07) 


(0.23 - 82.1) 


(0 - 100.5) 




shy/+BBy/ 





0.01 


0.11 


0.66 


5.09 






0-0.01 


- 0.04 


0.01 - 0.66 


0.07 - 3.02 


0.19-23.68 



Table 5: Synthesis of the numerical results for 3 objectives. Execution times are indicated in seconds. Symbol "-" means that the average execution 
time exceeds 30 minutes. 





n 


5 


10 


15 


20 


25 




Shy/ 





0.01 


0.02 






a = 0.4 




(0.9 - 2.6) 


(0 - 2.67) 


(0 - 4.33) 








shy/+BBy/ 





0.21 


10.31 










- 0.02 


0.02 - 0.57 


0.21 - 82.37 








Shy- 





0.01 


0.02 


0.05 




a = 0.5 




(0.63 - 2.03) 


(0 - 4.53) 


(0 - 4.9) 


(0 - 0.47) 






Shy/ +BBy/ 





0.12 


6.03 


669.09 








- 0.02 


0.01 - 0.72 


0.22 - 63.63 


10.84- 8779.1 






Shy/ 





0.01 


0.02 


0.05 


0.09 


a = 0.6 




(1.17-2.33) 


(0.07 - 7.8) 


(0 - 8.07) 


(0-3.1) 


(0 - 2.7) 




shy/+BBy/ 





0.11 


1.66 


84.87 


886.8 






- 0.03 


0.01 -0.32 


0.16-4.65 


6.15 -760.5 


12.46-5198.1 



Table 6: Synthesis of the numerical results for 4 objectives. Execution times are indicated in seconds. Symbol "-" means that the average execution 
time exceeds 30 minutes. 



tackled instances are more modest than for strictly decreasin g we ights. The only available alternative method should 
be however to resort to an exhaustive enumeration procedure 11611 . that would become intractable for 12 vertices. The 
results presented here are therefore a first step towards optimizing non-convex aggregation functions in multiobjective 
spanning tree problems. They make it possible to handle instances the size of which grows up to 25 vertices. 



6. Conclusion 

In this paper we have proposed several methods to solve the OWA-optimal spanning tree problem. One is based 
on a MIP formulation and is valid only when the weights of the OWA operator are strictly decreasing. We have 
shown that the use of a preprocessing phase, as well as a shaving procedure, makes it possible to considerably reduce 
the size of the problem, and therefore speed up the resolution. The numerical results prove the efficiency of the 
global resolution procedure (one is able to solve multiobjective instances with up to 100 vertices in reasonable times). 
Two branch and bound algorithms have also been presented, that work whatever weights are used. Despite a greater 
resolution time in the strictly decreasing case, they make it possible to tackle the non-monotonic case. This is a 
very challenging task, and the numerical results presented here are a first step in this direction. For future works, 
it would be worth further investigating optimization of OWA in this case, as well as some interesting variations of 



OWA, namely the non-monotonic OWA operator [42] (where negative weights are allowed) and the weighted OWA 
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n 


5 


10 


15 


20 


a 


= 0.4 


shy- 




(0.83 - 2.4) 


0.01 
(0.03 - 2.27) 


0.02 
(0- 1.17) 


- 






shy/ +BBy 


U 

- 0.02 


n a ^ 

0.01 - 2.98 


/in -7 a 

1.4-184.8 




a 


= 0.5 


shy 




(0.87- 1.93) 


0.01 
(0-1.6) 


0.02 
(0-0.1) 


0.06 
(0-0) 






Shy/+BBy, 


0.01 
- 0.02 


0.72 
0.05 - 4.45 


47.93 
1.27 - 528.6 


1674.9 
139- 8378.1 


a 


= 0.6 


Shy/ 




(0.47 - 1.2) 


0.01 
(0-1.47) 


0.04 
(0-1.4) 


0.05 
(0-0) 






shy+BBy 


0.01 
- 0.03 


0.87 
0.1 - 3.96 


40.12 
0.2 - 267 


760.9 
29.61 - 3611.1 



Table 7: Synthesis of the numerical results for 5 objectives. Execution times are indicated in seconds. Symbol "-" means that the average execution 
time exceeds 30 minutes. 



operator 13811 (where importance weights specific to each objective are allowed in addition to the weights of the 
OWA operator). Other promising research direction would be to propose a MIP formulation for a broader subclass of 
Choquet integrals (enabling to take into account positive and negative interactions between objectives), as for instance 
the class of fc-additive Choquet integrals 1 1311 . These research tracks are especially important because the search for 
a single best compromise solution is nearly the only operational approach when there are more than three objectives 
and the problem does not fit into the dynamic programming framework. 
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